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In a model with dynamical gluon mass, we investigate the 
bound states of two and three gluons via a Schrodinger equa- 
tion. The short distance potential is approximated by one- 
gluon-exchange while the long distance part is assumed to be 
of a breakable string. We estimate the masses and in particu- 
lar the sizes of low-lying bound states with no orbital angular 
momentum. By considering quantum-mechanical smearing 
of the gluon fields and normalizing to lattice results on M ++ 
and M 2 ++ , we find that the ++ glueball is rather small in 
size compared with the others. The fitted gluon mass is of 
order 600 to 700 MeV, which is reasonable. The 0" + , 1~~ 
and 3~~ three gluon glueball states are nearly degenerate, 
and their mass ratio with 2 ++ is largely independent of all 
parameters and consistent with lattice calculations. We esti- 
mate the mass of the 1 glueball to be around 3.1 — 3.7 GeV, 
which is close to the mass of J/tp and ip'. 

PACS numbers: 



I. INTRODUCTION 

During the last 20 years there has been much effort in 
trying to obtain a nonperturbative form for the gluon 
propagator. Perhaps one of the most interesting re- 
sult is that the gluon may have a dynamically generated 
mass jD. The existence of a mass scale, or the absence of 
a pole at k 2 = 0, is natural if one assumes that gluons do 
not propagate to infinity, i.e. these propagators describe 
confined gluons. The concept of massive gluon has been 
widely used in independent field theoretic studies, and 
examples about the consequences of massive gluons can 
be found in the literature The infrared behavior 

of the gluon propagator has also been studied numeri- 
cally jgj9|, and recent lattice computations give strong 
evidence for an infrared finite gluon propagator ]To| . 

The gluon self-coupling in quantum chromodynamics 
(QCD) implies the existence of bound states of gauge 
fields known as glueballs. Numerous technical difficulties 
have so far hampered our understanding of their prop- 
erties in experiments, largely because glueball states can 
mix strongly with nearby qq resonances. However, re- 
cent experimental and lattice studies of ++ , 2 ++ and 
h glueballs seem to be converging. All simulations 
agree that the lightest scalar glueball mass is in the range 
of 1500 — 1750 MeV, while the tensor and pseudoscalar 
masses are in the range of 2000 - 2400 MeV 0-0. It 
has been suggested |14[ that improved action lattice pre- 



dictions |nj agree well with the mass ratios of /o(1500), 
r/(2190), and /2(1980), which have exotic features that 
make them natural candidates for glueballs. All these 
states are seen in pp annihilation, central production in 
pp collisions, or J/ip — > 7 + X transitions |i"5|-|i"7|] . 

In this paper, we reopen the case of the potential model 
with massive constituent gluons, namely, the model of 
Cornwall and Soni |l8|,0. It is not our purpose to pur- 
sue the detailed theoretical bearings for gluon mass. But 
rather, we wish to explore how recent results in this field 
affect the potential model, and how the potential model 
can provide more insight on glueball properties. While 
the potential model has its limitations, it gives bound 
state solutions that have the advantage of providing in- 
formation such as the size of glueballs, an aspect which 
is rarely [^0| mentioned in the literature. In Sect. 2, we 
give the details of model description for low-lying bound 
states of 2-gluon and 3-gluon glueballs. We then use the 
variational method to estimate their masses and sizes in 
Sect. 3, where some smearing of gluon field is developed 
for the ++ case. Finally, we analyze our results and 
make some conclusions in Sect. 4. 



II. MODEL DESCRIPTION 

Although the dynamically generated mass m 2 (q 2 ) 
should be scale dependent, phenomenologically we shall 
treat m 2 (q 2 ) as constant for simplicity. The Lagrangian 
for the massive vector fields Af, is 
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where m is the effective gluon mass defined by Corn- 
wall [0 and others, and 



F^ = d^-d v A^+gf abc AlAl. 
The propagator for A^ is 
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In the nonrelativisic limit, we expand the massive gluon 
momentum and polarization vector to order |k| 2 , 
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gg scattering. 




where S = S1+S2 is the total spin of the 2-gluon glueball. 
As stated, we have transferee! the spin content of gluons 
to the wavefunction, and hence view iM fi as an operator 
acting on spin space. Note also that Si acts on ei and 
e3 while S2 acts on e 2 and e^. 

Using Eq. (0) , we arrive at the short distance potential 



Let us consider the system of a 2-gluon glueball. There 
are four diagrams contributing at tree level in the gluon- 
gluon interaction shown as in Fig. 1, corresponding to s-, 
t-, it-channel gluon exchange and the seagull. In the non- 
rclativistic limit, it turns out that all matrix elements are 
zcroth order in momentum except the contribution from 
s-channel, which can be removed both because of color 
and its second order nature in momentum. The short 
distance potential can be extracted from the tree-level 
Feynman amplitude of Fig. 1, 
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where q is the momentum transfer of the system. 

When extracting the potential, the "exchange" or 
"symmetric" diagrams are automatically taken care of 
by the properly symmetrized wave function for identical 
particle systems. Hence the relavent contributions are 



(13) 
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where sd stands for "short distance." 

The gluon-gluon interaction potential V s d in Eq. ( p^ ) 
is the Fourier transform of the tree level second order 
scattering diagrams of Fig. 1, but it cannot account for 
gluon confinement since it is of short distance nature. We 
must add a term to take into account such long distance 
effects. We add to V s d a string potential V^tr which is 
assumed to be spin independent, 
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(41 J p |2), (7) where (3 is related to the adjoint string tension Ka via 
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We define the spin operator S = S 2 , S 3 ) as 
(S ,fe )y = ~ie ljk , which satisfy 



In the potential V s t T , the color screening of gluons is 
brought about by a breakable string, that is, the adjoint 
string breaks when sufficient energy has been stored in 
it to materialize a gluon pair. This form of the string 
potential simulates the intergluonic potential as seen in 
lattice calculations 

We thus get the gluon-gluon potential relevant to two- 
gluon glueballs 
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i.e. S 1 , S 2 and S 3 are SU(2) group generators as desired. 
The spin operator can then be extracted via the relation 
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and placing the polarization vectors into the wavefunc- 
tion. We calculate the matrix elements in the center of 
mass frame, namely, Pi = pi = — P2, P/ = P3 = — P4, 
and the momentum transfer q = P3 — pi = Pf — Pi- 
After some simplifications, we obtain 
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where A is defined as 
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and is related to the strong coupling strength of the pro- 
cess. Note that the /-type constant in the last term of 
Eq. ([l3]) does not contribute when contracted with the 
normalized color wavefunction 
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We are left with three parameters: effective gluon mass 
m, string breaking parameter and adjoint strong cou- 
pling constant A. We take M the conservative range 
j3 ~ (1 ± 0.7), while A is determined by |1 
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for SU(iV) group. For TV = 3, taking m ~ 600 MeV, A ~ 
350 MeV, one gets A ~ 1.4. For A varying from 250-400 
MeV, A varies in the range 1.1 — 1.6. 

For the case of 3-gluon glueballs, we assume that the 
constituent gluons interact pair- wise (HI. Thus 
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We note that the contribution from s-channel gluon ex- 
change can still be ignored because of its second order na- 
ture in momentum. The factor of one-half for the string 
potential is because one needs to pull three (and not six) 
gluons from the vacuum to screen the three gluons that 
are originally there in the glueball. 

For the low-lying bound states with relative angular 
momentum ijj = for each pair of gluons in the 3-gluon 
system, the normalized color singlet wavefunction is ei- 
ther symmetric or antisymmetric: 
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After contracting these color wavefunctions with the 
short distance potential of Eq. (|l^) , we get the pair- wise 
potential for 3-gluon glueball case 
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for d- and /-type, respectively. 



III. GLUEBALL MASSES AND SIZES 

In this paper, we consider only the case of L = 0, 
hence ignore the spin-orbit and tensor terms in Eqs. ( |l6| ) 
and (g3|). The 2-gluon and 3-gluon glueball systems are 
discussed separately. 



A. Two-gluon Glueballs 

For 2-gluon glueballs with L — 0, we have only J PC = 
++ ,2 ++ states. The Hamiltonian is 
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We immediately notice one serious problem: when S = 0, 
one has an attractive (^-function term and the Hamilto- 
nian is unbounded from below. This "maximum attrac- 
tion channel" in ++ could be related to the gluon con- 
densation that triggers confinement. 

In Ref . |ll| , the (5-function term was treated as a per- 
turbation. For our study, we propose a physical solution 
by smearing the gluon fields, that is, we replace the 5- 
function by the smearing function 
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which approaches <5 3 (r) for k — > oo. Using the varia- 
tional method with trial wavefunction ip{ r ) ^ e ~ a m r > 
we illustrate in Fig. 2 the smearing dependence of 2- 
gluon glueball masses for A = 2 and j3 = 0.3. 

We see from Fig. 2 that the mass of the 2 ++ glueball 
converges rapidly to 3.2 m for k > 1, illustrating good be- 
havior since the ^-function term is repulsive, but the ++ 
mass decreases monotonically until it becomes negative 
for k ^ 3.8. To illustrate what is happening when the 
attractive S becomes operative for larger k (less smear- 
ing), we plot in Fig. 3 the root-mean-squared radius 
r rms = yj (r 2 ) of ++ and 2 ++ glueballs as a function of 
k for the same values of A and (3. It is clear that, while 
the size of 2 ++ glueball stabilizes for k > 1, the radius for 
++ glueball drops monotonically with its mass, which in 
turn drops monotonically with increase of k. 
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FIG. 2. M 0+ + (solid) and M 2 ++ (dashed) vs. smearing 
parameter k for X—2, /3—0.3. In large k limit, one recovers 
the 5-function potential which drives M ++ negative.. 
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FIG. 3. Glueball radius for 0++ (solid) and 2++ (dashed) 
glueballs vs. smearing parameter k. 

One may be tempted to use Fig. 2 and argue that, 
since for 1 < k < 3, both 2 ++ and ++ glueball masses 
are relatively stable, hence they have approximately of 
the same mass. However, this is not what is observed 
on the lattice nor as suggested by experiment. Since we 
cannot claim to know how to determine the value of k in 
Eq. (p6|), we use the converging experimental and lattice 
results to fit for k. We take M ++ = 1730 MeV, and 
M 2 ++ = 2400 MeV [OJT1, hence the mass ratio 
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= 1.39. 



o++ 



(27) 



Using Eq. ( p7| ) to determine k, we find that it depends 
mainly on A, but is almost independent of (3. This is to 
be expected since the need for smearing comes from the 
short distance potential. With the mass ratio fixed, we 
find a smaller A (and a greater (3 in general) can accom- 
modate a larger k value, which is reasonable. The typical 
value of k falls in the range of 2.3 to 4.3, which ensures 
that the 2 ++ mass is stable. 
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FIG. 4. M 0++ and M 2++ vs. (3 for A = 1.5, 2.0 and 2.5 as 
indicated by the symbols "+" , "o" and "x" , respectively. To 
guide the eye, the points are linked by solid and dashed lines. 
The mass ratio is held fixed by Eq. fej) 
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FIG. 5. Glueball r rms radius for ++ and 2 ++ glueballs vs. 
/3 for A = 1.5, 2.0 and 2.5 as indicated by the symbols "+" , "o" 
and "x", respectively, and linked by solid and dashed lines. 

With k determined for given A and (3, the mass and 
size of 2-gluon glueballs can be calculated. Since M ++ 
is lighter than M 2 ++ from experiment/lattice input, the 
++ glueball is rather small in size compared to the 2 ++ 
glueball. This is reasonable because of the attractive (re- 
pulsive) 8 function potential for ++ (2 ++ ). Except for 
this smallness in size of ++ glueball, which is dominated 
by A part of potential, the ++ mass, and the mass and 
size of the 2 ++ glueball are all more dependent on f3, the 
string part of potential. We plot the masses (in units of 
m) and sizes (in units of 1/m) vs. (3 for three different 
values of A in Figs. 4 and 5, respectively. These figures 
illustrate the range of uncertainties within the model. 
We find that the masses increase with increasing (3 or 
decreasing A; however, increasing j3 or decreasing A will 
increase the size of 2 ++ glueball but decrease the size of 
++ glueball. We will give a more detailed discussion on 
this later. 

As mentioned already, the 2 ++ glueball mass is stable 
and almost independent of k for k > 2, which holds for 
all solvable (A, (3) parameter space. Hence, taking the 
2 ++ glueball mass range of 3.4 m to 3.9 m from Fig. 4 
and the lattice result of M 2 ++=2400 MeV, we estimate 
the effective gluon mass to be 0.6 — 0.7 GeV, which agrees 
well with the gluon mass ~ 0.66 GeV needed |,| to ex- 
plain the photon spectrum in radiative J/ip decay. With 
gluon mass m determined, the typical size of 2 ++ and 
++ glueballs can be read off from Fig. 5, which are in 
the ranges of 0.8 — 1.1 fm and 0.1 — 0.2 fm, respectively, 
which is in good agreement with the results obtained 
on a lattice using the source method. For a more direct 
calculation of the ++ glueball mass on the lattice, one 
would need relatively fine lattice spacings pU. It would 
be interesting to see if our result of small 0" size could 
be further replicated on the lattice. 



B. Three-gluon Glueballs 



For 3-gluon glueballs, we introduce the center-of-mass 
and relative coordinates: R = (ri + r 2 + r 3 )/v3, ri 2 = 



(ri — r 2 )/v / 2 ) an d r = \/2/3(i"i + r 2 — 2r 3 ). Since we 
consider only pure s-states, one has the pair potential 
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where S pa i r is the spin of any pair of gluons in a 3-gluon 
glueball system. The Hamiltonian for this system is then 

1 



Hd.f = 3m - J-V? - — VJ ia 
2m 2m 12 

+ Vi,/(ri2) + Vd,/(r 23 ) + ^,/(r 3 i). 



(29) 



Since the glueball wavefunction must be symmetric 
with respect to ri, r 2 , and r 3 , the contributions of the 
three pair potentials are the same. Hence the spin and 
spatial parts of the glueball wavefunction are indepen- 
dent and 



J2 S^air = (Si + S 2 ) 2 + (S 2 + S 3 ) 2 + (S3 + Si) 2 
= (S 2 + S 2 + S 2 ) + (S 1 + S 2 + S 3 ) 2 . 

The Hamiltonian above can then be simplified as 
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where Stotai is now the total spin of the system, i.e. J. In 
Ref . JTsfl , one introduced an additional quantum number 
called S'pair in J PC (S pa i T ) which is clearly not adequate, 
resulting in a spurious 1~~ state. 
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FIG. 6. Three-gluon glueball masses M -+, M x ~- and 
M 3 -- vs. /3, as indicated by the symbols "+", "o" and "x", 
respectively. 
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FIG. 7. Glueball radius vs. (3 for 0"+, l -- and 3~ 
3-gluon glueballs, as denoted by the symbols "+", "o" and 
"x", respectively. To guide the eye, the points are linked by 
solid lines. 

There are three possible pure s-states, with quantum 
numbers J PC = h , 1 , and 3~~. For these glueballs, 
all the 8 terms are repulsive, hence their masses are au- 
tomatically bounded from below, and there is no need to 
smear the 8 function. We estimate the masses and sizes 
of these glueballs using the variational method with the 
wavefunction ip(r) oc e~ a m ( r + r i2) ; which is symmetric 
with respect to ri, r 2 , and r3. We hnd that the masses 
and sizes are almost A independent, hence we plot them 
only vs. (3 in Figs. 6 and 7, respectively. Note that for 
the ^ quantum number there is also a 2-gluon glueball 
in L = 1 state, which is expected to have mass lower 
than the 3-gluon state discussed here fL8| ]. Mixing be- 
tween the 2-gluon and 3-gluon states should lead to level 
repulsion and raise the 3-gluon state above the two 
other states. 

From Fig. 6, we find that the masses of the three lowest 
lying 3-gluon glueballs are within 0.1 m (or 100 MeV) of 
each other, which holds for any value. In other words, 
they are nearly degenerate, which we will discuss further 
in the next section. We note that the ratio with 2 ++ 
mass is around 1.53, almost independent of f3 and A, and 
agrees with naive constituent counting. Scaling from the 
lattice result of M 2 ++ = 2000-2400 MeV the mass range 
of these glueballs is 3.1 — 3.7 GeV, right in the ballpark 
of J/tp and ip' masses. Their sizes are only slightly larger 
than the 2 ++ 2-gluon glueball, and fall into the range of 
0.9 - 1.3 fin. 



IV. DISCUSSION AND CONCLUSION 

We approximated the attractive ^-function by a smear- 
ing function D(r) in Eq. (|26|). This is physically rea- 
sonable since the gluons in a glueball can not have very 
large relative momentum that corresponds to very short 
distance separation. Through the variational method, 
besides the constituent terms 2m and 3m, one can check 
four sources for the glueball mass: kinetic energy, Yukawa 
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and delta function potential, and string energy. It is clear 
that smaller glueball size implies a greater kinetic energy 
contribution to its mass. The Yukawa term is always at- 
tractive and contributes negatively to the glueball mass. 
It is proportional to the adjoint strong coupling constant 
A, Eq. dl7| ), and becomes more negative for smaller glue- 
ball size. The delta-function term, replaced by the smear- 
ing function D(r) of Eq. (|26|), is also proportional to A, 
and becomes stronger for smaller glueball size. As for 
string energy term, a bigger glueball can store more en- 
ergy hence a larger glueball mass. A greater string ten- 
sion (3 also stores more energy in the string. Knowing 
about these four sources, we can see how the mass and 
size of a glueball comes about. 

Let us first consider masses. For the 2 ++ glueball, 
the kinetic energy is relatively small due to its rather 
large size compared to the ++ glueball. On the other 
hand, there is a cancellation between A-dependent attrac- 
tive Yukawa potential and repulsive delta- function terms. 
The resulting value is always negative, and will further 
cancel against the kinetic energy term. Hence the main 
contribution to the mass comes from the string energy, 
and the stronger the string tension (3, the heavier the 
mass. For fixed string tension, a larger A value gives a 
stronger cancellation between the remnant of short dis- 
tance potential and the kinetic energy, resulting in a 
smaller contribution to the glueball mass. In other words, 
the 2 ++ glueball mass increases with increasing (3 or de- 
creasing A, as can be seen in Fig. 4. 

For the ++ glueball, the kinetic enengy term is rel- 
atively large because of the small size, which in turn is 
brought about by the attractive 8 function, and there is 
strong cancellation between these two terms. Since we 
have fixed the mass ratio with 2 ++ , it turns out that the 
kinetic energy overcomes the attractive 8 function, and 
subsequently cancels against the attractive Yukawa term. 
The net result could be positive or negative depending 
on smaller or larger A value. The main contribution still 
comes from the string term but it cannot be as strong as 
in the case of 2 ++ because of the small size. We stress 
again that in Fig. 4 we have held the 2 ++ and ++ ratios 
fixed according to Eq. ( ]27j ) . 

We turn to the consideration of sizes, i.e. understand- 
ing Fig. 5. We find that the glueball size increases with 
mass for 2 ++ glueball. This is easily understood in terms 
of a repulsive 8 and the dominance of string energy. For 
0++ glueball, its size is sensitive to A, which can be un- 
derstood as coming indirectly from smearing. As a con- 
sequence of imposing Eq. (^7j), the value for k depends 
mainly on A and can become larger for lower A. For larger 
k, the smearing function D(r) approaches the delta func- 
tion and is more attractive. Thus, paradoxially, a smaller 
A ("less attractive") is more able to "pull in" the glue- 
ball, leading to decrease in size. However, the size of the 
++ glueball is almost (3 independent precisely because of 
its small size, hence insensitive to the string energy. We 
note that the smallness of the ++ glueball stretches the 
applicability of our relativistic expansion. It is therefore 



amusing that Figs. 4 and 5 are in rather good agreement 
with the findings of Ref. [^(J . 

For three-gluon glueballs, there are three possible pure 
S states, namely J PC =Q 1", and 3 . For these 
glueballs, all the 8 terms in the potential are repulsive, 
hence no smearing is needed. There is a A-dependent can- 
cellation between the Yukawa and delta-function terms. 
Interestingly, both terms become stronger for larger to- 
tal spin (or J), and the cancellation hides the spin effect, 
resulting in the masses being nearly degenerate. In other 
words, their masses and sizes are almost A-independent 
and depend basically on (3, just like 2 ++ case. It is rather 
intriguing that the mass difference of the three glueballs 
are within 100 MeV of each other. But as we have men- 
tioned earlier, the ^ state would become heavier via 
mixing with the L — 1 two-gluon state. 

Comparing Figs. 4 and 6, we note that the 3-gluon 
glueball masses are about 1.5 times larger than 2 ++ , 
largely independent of (3 and A. Taking M 2 ++ between 
= 2267 ± 104 MeV |§ and 2400 ± 25 ± 120 MeV [[n), 
one finds M x -- in the range of 3.5 — 3.7 GeV, close to the 
ip' mass of 3686 MeV. We note that the proximity of the 
1 glueball to if) 1 may be called for from comparsion of 
J/i[i and ip' two body hadronic decays Jp2j. However, if 
the / 2 (1980) state is the 2++ glueball ]u|, tnen wc nnd 
Mi-- is of order 3 GeV or closer to the J/tp, where an 
older proposal [^3) of l -_ glueball (called O) could be 
behind the rather sizable strength of J/ip — > pw ~ 1%. 
At the moment, improved action lattice results jll|] find 
1~ _ glueball masses heavier than those discussed here. 
We urge further refined, dedicated studies to help clarify 
the phenomenology. Direct search methods for the lowest 
lying 3-gluon glueballs were discussed in Ref. (T^| , which 
should also be brought up to date. 

In conclusion, we investigate the bound states of two 
and three massive gluons with a Schrodinger equation. 
We calculate the short distance potential from one gluon 
exchange, and give arguments for the long distance con- 
fining potential. We calculate glueball masses and sizes 
using variational method. By considering the effect of 
smearing of the gluon fields, we find that the size of 
++ glueball could be rather small compared with oth- 
ers. The other glueball masses are stable with respect to 
such smearing. Using the converging experimental and 
lattice results for ++ and 2 ++ glueballs, we estimate the 
effective gluon mass to be 0.6 to 0.7 GeV, in agreement 
with phenomenological results. The typical size of 2 ++ 
and ++ are of order 1 fm and 0.1 to 0.2 fm, respectively. 
This means that to extract the ++ glueball size on a 
lattice, one would need rather fine lattice spacings. It 
would be of interest to see if our result would be repli- 
cated on the lattice. For 3-gluon glueballs, their sizes are 
also estimated to be in the range of 0.9-1.3 fm, similar 
to that of the 2 ++ state. We find that the three lowest 
lying 3-gluon glueballs to be largely degenerate. In par- 
ticular, the mass of the 1~~ glueball is in the range of 
3 to 3.7 GeV, which is consistent with arguments from 
phenomenological point of view. 
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